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We have investigated the role of proton-neutron pairing in 
the context of the Quasiparticle Random Phase approxima- 
tion formalism. This way the neutrinoless double beta decay 
matrix elements of the experimentally interesting A= 48, 76, 
82, 96, 100, 116, 128, 130 and 136 systems have been cal- 
culated. We have found that the inclusion of proton-neutron 
pairing influences the neutrinoless double beta decay rates sig- 
nificantly, in all cases allowing for larger values of the expec- 
tation value of light neutrino masses. Using the best presently 
available experimental limits on the half life-time of neutrino- 
less double beta decay we have extracted the limits on lepton 
number violating parameters. 
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I. INTRODUCTION 

Among the exotic processes the neutrinoless double j3- 
decay (0z//3/3-decay) 



(A, Z) ->(A,Z± 2) 



(1) 



has been sought experimentally for about half a century 
but it has not yet been observed Q . Its observation will 
undoubtedly be a signal of interesting physics beyond the 
standard model of electroweak interactions. First of all it 
will demonstrate the breakdown of lepton number conser- 
vation which, being associated with a global, not gauge, 
symmetry is expected to be broken at some level. It will 
also give us useful information about the neutrino mass if 
it is in the region < 1 eV. Finally, and most importantly, 
it is the best process, if not the only one, to settle the 
issue of whether the neutrino mass eigenstates are of the 
Majorana type, i.e., whether the particle coincides with 
its own antiparticle (a la ir°), or of the Dirac type (a la 
K°). 

It is, therefore, not surprising that the experimental 
searches M- [n5| for the above process have persisted with 
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great devotion up to the present day and have lead to the 
unbelievably long life time limit Q 

T°J 2 > 5.6 x 10 2 V 

From this limit, in conjunction with calculations of the 
nuclear matrix elements involved, the limit | < m„ > | < 
0.8 eV has been extracted for the average light neutrino 
mass. It is interesting to remark that neutrino masses in 
this neighborhood can constitute candidates for the Hot 
Dark Matter (HDM) component required for the under- 
standing of the large scale structure of the universe as 
indicated by the KOBE data 

The analysis of the OV/3/3-decay data, if and when they 
become available, is unfortunately not going to be sim- 
ple. Depending upon the extension of the standard model 
assumed there are many mechanisms which can lead to 
pro cess (^) which can interfere with one another (see refs. 
Hl8"f - [62) for reviews). In those mechanisms in which the 
exchanged particle is light (e.g., light Majorana neutrino) 
the effective transition operator is of long range (essen- 
tially coulombic) with or without spin dependence. If 
on the other hand the exchanged particle is heavy (e.g., 
heavy Majorana neutrino, exotic Higgs Scalars etc.) the 
effective operator is of a short range and in the presence 
of short range correlations, one must be careful not to 
ignore (ll| the finite size of the nucleon (~ 0.8 fm). In 
this case the size of the operator is set by the size of the 
nucleon. It is obvious that the two cases involve different 
nuclear physics and as we shall see later they lead to the 
extraction of different parameters of the gauge theory. 

It is clear from the above discussion that all nuclear 
matrix elements must be computed reliably. This is eas- 
ier said than done, however, since the nuclear sytems 
which can undergo double /3-decay, with the possible ex- 
ception of the A = 48 system, are far away from closed 
shells and as a result have complicated structure p8|- 
pl[ . In the shell model approach it is clearly impos- 
sible to construct all the needed intermediate states of 
the nucleus which is one unit of change away. Since, 
however, the energy denominators are dominated by the 
virtual neutrino momentum, rather than the nuclear ex- 
citation energy, the construction of all these states can 
be avoided using a version of the closure approximation 
|i~8f - But even then it is quite hard to construct the 
wave functions of the initial and final nuclei employing 
a full shell model basis. Thus a weak coupling scheme 
has been employed, starting from products of neutron 
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and proton wave functions and employing truncations 
according to the energies of the unperturbed proton and 
neutron states. In any case it is quite clear that it is very 
hard to substantially improve these old calculations 

It is thus not surprising that other approximation 
schemes have been employed. The most prominent 
among them has been the Quasiparticle Random Phase 
Approximation (^6|- |}4| (QRPA). In this scheme the 
construction of the intermediate states is unavoidable, 
but not extremely hard. In the first step one constructs 
the intermediate states \mJM > as proton-particle neu- 
tron - hole excitations built on the ground state of initial 
nucleus. In the second step one views the intermediate 
states |m JM > as neutron-particle proton-hole, built on 
the ground state of the final nucleus. These intermediate 
states are expressed as two quasiparticle states and one 
makes proper adjustments for the fact that the two sets 
of states \mJ m > and |m J m > are not orthogonal. 

In the QRPA approximation there is no need to in- 
voke the closure approximation p(J- ps| . In fact it was 
possible to use QRPA to explicitly check how well the 
closure approximation works. It was rewarding to find 
that it works quite well except in those situations when 
the matrix elements are unusually suppressed. 

As we have previously mentioned a number of QRPA 
calculations (24|- |58| for almost all nuclei of practical 
interest in O^/3/3-decay have been performed. In all 
such calculations the two quasiparticle states were of 
the proton-neutron variety. Proton-neutron (p-n) pair- 
ing correlations |39|- |TJ| had been neglected. Such cor- 
relations were recently found, however, to be important 
p4| fUj[ in the evaluation of the nuclear matrix elements 
entering 2i//3/?-decay 



(A, Z) -> (A, Z + 2) + e~ +e~+P e + , 



(A, Z)^{A,Z-2) 



(2) 



(3) 



which proceeds only via the 1 + intermediate nuclear 
states. 

It is the purpose of the present paper to investigate 
what effect, if any, the p-n pairing has on the Oi^/3/3-decay 
matrix elements. To this end we will repeat and extend 
our previous calculations [p6|-p8t to cover most of the nu- 
clear targets of experimental interest ( 48 Ca, 76 Ge, 82 Se, 
96 Zr, 100 Mo, 116 Cd, 128 Te, 130 Te, 136 Xe). 



II. LEPTON VIOLATING PARAMETERS AND 
ASSOCIATED NUCLEAR MATRIX ELEMENTS 

As we have mentioned in the introduction there exist 
many mechanisms which can lead to 0v(3(3 decay some 
of which are exotic (exotic Higgs scalars, supersymmet- 
ric R-parity violating interactions etc. ]lql ). The most 
popular scenario is one which involves intermediate mas- 
sive Majorana neutrinos. We will concentrate on the last 



mechanism in this work even though some of the nuclear 
matrix elements computed may be used for the more ex- 
otic mechanisms as well. 

We will distinguish the following two cases: 
i) Both leptonic currents are of the same chirality, i.e., 
both left handed or both right handed. Then out of the 
neutrino propagator one picks the part 



(4) 



,2 ' 



which for light neutrinos yields an amplitude propor- 
tional to 



< m v >- 



Em 



(ii) 

ej 



(5) 



exp(iaj) is the CP eigenvalue of the neutrino mass eigen- 
state |z/j >. For light neutrinos the contribution of the 
right handed currents is negligible. For heavy neutrinos 
the amplitude becomes proportional to 
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(7) 



for left and right handed currents respectively, exp(i^j) 
and Mj are the CP eigenvalue and the mass of the heavy 
neutrino mass eigenstate \N\ >. The sub-matrices U 1 - 11 ^, 
[/(i2) 5 JJ^ 21 \ £/( 22 ) are the parts of the unitary matrix 
which connect the weak eigenstates v^ c ), a — e, fi, r, 
with the mass eigenstates (ui,N{), i = 1, 2, 3 (see e.g., 
rcf. JL8|). The upper index "1" stands for the light and 
"2" for the heavy neutrinos of left-right symmetric mod- 
els. For each case one encounters a spin independent 
operator (vector x vector or Fermi) and one which is a 
spin dependent scalar (axial x axial or Gamow- Teller), 
yielding the matrix elements Mp and A/qt respectively. 

ii) leptonic currents of opposite chirality. 
Then the relevant part of the intermediate neutrino prop- 
agator is 



1\ 



,2 ' 



(8) 



Since the right-handed currents are expected to be sup- 
pressed, this contribution is expected to be significant for 
light neutrinos. Thus the amplitude does not explicitly 
depend on the neutrino mass. One can now extract from 
the data two dimensionless parameters Q 



A 



(9) 



(10) 



where 



2 



K = (MJM R ) 2 , 
- = tg £ (mixing). 



(11) 
(12) 



Ml and Mr are respectively the masses of the vector 
bosons W L and iy R associated with the left and right 
handed interactions. £ is the W h — W R mixing angle. 
Vul is S iven b y @ 



(13) 



with 



^ = ^2T + (i)T + {j)u(i,j)g m (qi,qj), 



3m(qi,qj) 



47ri?o 1 <5(qi + qj) 
(2tt)3 2 v / 2A m (g i V2)' 



(16) 



(17) 



R = r a A 1/z (r ~ 1.1. fm) (nuclear radius), (18) 



The A-term arises when the chiralities of the hadronic 
currents match those of the leptonic currents, i.e., they 
are of the Jl — Jr combination. The ?7-term arises when 
the two hadronic currents are of the same chirality, i.e 
of the Jl — Jl type (this is possible due to the W-boson 
mixing). For the extraction of A one should know five 
matrix elements i.e., (Mp w , MstoO which arise from 
the time component of the propagator of Eq. (||) and 
(Mp, Mq T , M t (tensor)) arising from the space compo- 
nent of that propagator. Due to the different energy de- 
pendence Mpui , Mgtw , Mp and M' GT are different from 
the two matrix elements Mp and A/qt encountered in 
the light neutrino mass mechanism. We will see, how- 
ever, that to a good approximation 



M F 



M 



GT 



Mp. 



Mq t 



M, 



GTw 



(14) 
(15) 



The situation is a bit more complicated in the case of 
the ?7-term since one encounters two additional matrix 
elements Mp and Mr. The first arises from an opera- 
tor which is antisymmetric both in the spin and angu- 
lar momentum indices jl8| |- pl| . The second matrix ele- 
ment arises from the momentum dependent terms of the 
hadronic current pEf (weak magnetism etc.). The con- 
tribution of the momentum dependent term is normally 
a small correction. This is not true in this case, however, 
since due to the structure of the propagator of Eq. (||), 
the standard leading non- vanishing term is proportional 
to the average lepton momenta jD . 

One must note, furthermore, that the kinematics are 
different, reflecting the difference between (||) and (^). 
Thus the coefficients entering the various combinations 
of nuclear matrix elements in the decay rate are energy 
dependent. As a result the relative importance of the 
various nuclear matrix elements may vary from nucleus 
to nucleus (depending on the available energy). 

In the case of heavy neutrino intermediate states one 
encounters two matrix elements Mhf and Mhgt which 
differ from the above matrix elements Mp and Mgt due 
to the fact that the radial part of the relevant operator 
is short-ranged. 

The nuclear matrix elements mentioned above are as- 
sociated with a set of transition operators which in mo- 
mentum space can be cast in the general form (see refs. 
|6| MM for details) 



A m (g) = vV + m l [ £ m + \J q 2 + ml 



(19) 



where m„ is the mass of the virtual neutrino and e m a 
suitable energy denominator [see Eq. ( fi6] ) below]. 

The most important matrix element Mqt is associated 
with a i ■ <?2, i.e., 



Similarly 



wgt(1, 2) = a x ■ <? 2 < — ► M G t- 



cj f (1,2) = 1 < — ► M F 



(20) 



(21) 



The matrix elements Mfh,Mgth associated with 
heavy neutrino are related to the operators wfh and 
wqth where 



wfh(1,2) 



w G th(1,2) 



m G m P 



M 



FH, 



m n m r 



-<J\ ■ U2 



Mo- 



th- 



Notice that for heavy neutrino 

ml 1 1 



m e m p A m (q) m e m p 
(independent of momentum) to be compared with 

2 (for light neutrino). 



A m (g) 



(22) 



(23) 



(24) 



(25) 



For processes which do not explicitly depend on the 
neutrino mass (j'l — Jr. interference) we encounter the 
operators cjfu and wgtu which differ from wp and wqt 
by the inclusion of the extra kinematical factor 5 m {q\\/2) 
with 



y/q 2 + ml 



+ \/q 2 + ml 



(26) 



One also encounters the matrix elements Mp and M GT 



associated with the operators 
Wp(l,2) = -2gi-V x 



M p , 



lo' gt (1, 2) = -1B X ■ a 2 qx ■ Vi < — > M GT , 



(27) 



(28) 
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and the matrix element Mp associated with 



and 



w T (l,2) = -<ti • a 2 q\ ■ Vi - 3\ ■ qia 2 ■ Vi 



Xi± = ±3xf + Xgt ~ 6 Xt> 



C?l • Vi<7 2 • ?! 



M T . 



(29) 



As it has been mentioned above in the case of j'l — Jr but 
Jl — Jl combination we encounter two additional matrix 
elements namely M P and Mr . M p is associated with the 
operator w P (l, 2) : 

w P (l,2) = - B 2 ).{l x -l 2 ) <—> M P , (30) 
while Mr is associated with the operator oj r : 

> Mr, (31) 



^ fly 2 M r 2 , ,r , ,Ri 



<7 A m e m p a 3 

with /i = fi p — // n -|_i = 4.7 (a is the oscillator length) and 
u§{l,2) = {5 1 -a 2 )q 2 a 2 , (32) 

2) = (ai ■ qi <? 2 ■ qi - ^<?i • a 2 )q 2 a 2 . (33) 

In the expression for the Oi^/3/3-decay life time vari- 
ous combinations of the above nuclear matrix elements 
appear. These will be given in units of M GT and be de- 
noted [^9) by x- First those associated with the mass 
mechanism 

<nv> )+< mp 

m e Mn 



with 



X K = (n 2 + e 2 ) < -f- > R Xh , 



Xh = (( — ?M FH - M GTH )/M GT . 
5a 



(35) 



(36) 



(37) 



Secondly those not connected with the mass mechanism: 
Xfw; Xgt^, Xr, Xi±, X2± Xf> Xgt> Xt> Xp where 



_ (£V\2 M Fu 

3a m gt 



M, 



XGTa 



XR 



GTu 



M GT ' 



Mqt 



(38) 
(39) 
(40) 



X2± = ±X 



XGTa 



9 Xl± 



(41) 



(42) 



in an obvious notation (x F = Mp/M G T etc.). In the limit 
in which the energy denominator e m can be neglected we 
obtain 



Xf = Xf = XFu 



XGT = XGT = XGTa 



1. 



(43) 



(44) 



As we have already mentioned the closure approxima- 
tion was not employed in our approximation. In writing, 
however, the expression for e m we made the standard 
approximation of replacing the electron energy by an av- 
erage value. Furthermore the lepton wave functions were 
taken out of the nuclear integral incorporating their effect 
via a distortion factor (see refs. [p^| , p^| ). 

In avoiding the closure approximation [p6|- p8p9| ] the 
momentum space representation was found extremely 
useful. As it was shown in ref. ]39| by exploiting the 
properties of the harmonic oscillator wave functions it 
was possible to express the energy dependent radial in- 
tegrals of each type of operator in addition to the energy 
in terms of only two parameters n and £. The parameter 
£ takes values (£ — 0,1,2 for scalar, vector and tensor 
rank respectively) . n also takes few values limited by the 
number N of oscillator quanta involved in the two par- 
ticle wave function of the interacting nucleons (n < 2N) 
regardless of the number of single particle configurations 
employed [j39| . In the systems we studied n < 10. 



III. THE QRPA FORMALISM WITH 
PROTON-NEUTRON PAIRING 

The QRPA formalism employed in the /3/3-decay, a 
modification of the usual RPA formalism that involves 
a change of charge by two units, is fairly well known 
p6[- p7| . So we will briefly mention its main features 
here. (3(3 decay is viewed as a two step process. In the 
first step a proton particle neutron hole is created on the 
initial state |0^ > and associated with the intermediate 
states \m,JM > of the (A, Z + 1) nuclear system. In 
the second step, rather than considering an additional 
proton particle - neutron hole acting on the excited in- 
termediate states, we consider a neutron particle-proton 
hole on the final state |0^ > leading also to states of the 
(A, Z + l) nucleus labeled \mJM >. The states \m.JM > 
and \mJM > are obtained by solving the well known 
QRPA equations. Unfortunately they are not orthogo- 
nal to each other. The nuclear physics dependence of 
/?/3-decay is essentially contained in the matrix element 
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< O f + || u k (a+ a nd ) || mJ > 

x < mJ\m.J >< mJ | u k (a p3u a nh ) || 0, + >, 



B^(a 7 b,JM) = (-1) J " M B (a,6, J,-M), (51) 



(45) 



where u\(a + a) are the usual tensors of rank k (k = 
J here) built of protons and neutrons as indicated. 
u\(a + a) = (—l) k ~ x u k _ x (a + a) is the time reversed op- 
erator of u k (a + a). a, b, c, d designate all the single parti- 
cle quantum numbers except for the angular momentum 
projection m a i.e., a {n a ,i a ,j a ) etc. The overlap is 
necessary since, as we have mentioned above, these inter- 
mediate states are not orthogonal to each other. Hence 
the energy denominator e m encountered in the previous 
section is given by the prescription 



E m — < E c > —E[ — E n 
\(E m + E w )-^(E i+ E { ) 



(46) 



Here, fi m and f2m are the QRPA energies of the excited 
states \mJM > and \rnJM > calculated from the ground 
state energy of the initial and final nucleus, respectively. 

The next step consists in expressing the tensor opera- 
tors u\(a + a) in terms of quasiparticles |26|| i.e., 



pam a 



n a m a 



_ ( u p {a) 
Up (a) 



-«p(a) 
■Up(a) 



"lam a 
^lam. 



u n (a) -v n (a) 
v n (a) ttn(a) 



2am a 

h 

J 2&m a 



(47) 



(48) 



where b + and b are the construction and destruction op- 
erators for quasiparticles. The tilde ~ indicates the time 
reversed states a a m a = (— )j»~™ a a a _ nia etc. Clearly since 
protons and neutrons do not mix the indices 1 and 2 
can be identified with protons and neutrons respectively. 
The parameters u and v are the occupation probabili- 
ties obtained by solving the standard BCS equations for 
the initial state. A similar set of equations for the final 
states (A, Z+2) yields the occupation probabilities u and 
v entering the matrix element in the left of Eq . (BaV 
The QRPA states \mJM > are of the form |§§~|l) 

\mJM> =QT M + \0t> RPA 

= J2 [X£(a,b,J)B+(a,b,J,M) 

a, b 

- Yg(a, 6, J)B 12 {a, b, J, M)]\0+ >, (49) 



where 



B »„(a, b , J, M ) = n(fJ,a, vb) 



E 

m a ,m b 



r JM h+ h+ (it v 

^j a m a j b m b >am a u l /bm VP, " 



1,2), 



(50) 



and 



n(fia, vb) 



(1 + 

(l+^(5ab) 3 / 2 ' 



(52) 



The forward- and backward- going amplitudes Xf^ and 
and the energies of the excited states fi m are ob- 
tained by solving the QRPA equation for the initial nu- 
cleus (A,Z) By performing the QRPA diag- 
onalization for the final nucleus (A, Z+2) we obtain the 
amplitudes X\ % , and eigenenergies Jig- of the QRPA 
state \mJM >. 
One can show that 



<mJ\\ M k (a+a nb ) || 0+ > 
= 4j V2J+T [u p (a)v n (b)X™(a,b, J) 
+ v p (a)u n (b)Y£(a,b,J)}, 



(53) 



< 0+ || u k (a+ c a nd ) || mJ > 
= <5 k j V2JTT [v p (c)u n (d)Xf 2 (c, d, J) 
+ u p (c)v n {d)Y™{c, d, J)}. 

The overlap integral takes the form 

< mJ\mJ >= 5}*S(a, b, J)Xf 2 {a, b, J) 

a, b 



(54) 



Y3(a,b,J)Y3(a,b,J)]- (55) 



Once proton-neutron correlations are introduced the 
quasiparticle labels can no longer be identified with pro- 
tons and neutrons, but the y b ecome mere labels. The 
transformation matrix [^0|- |f6| is generalized to the 4x4 
matrix 



pa m 
+ ' 
nam a 



^pam a 
y ^nam a 



( u lp (a) u 2p (a) -Vi p (a) 

ui n (a) u 2n (a) -«in (a) 

wip(a) v 2p (a) u lp (a) 

\ Wi n (a) V2n(a) ui n (a) 



-v 2p (a) 

-V2 n (a) 
U2n(a) 

h + 

u 2am a 
&lam a 
V &2am a / 



(56) 



independently of the angular momentum projection 
quantum number m a . 

The columns of the above matrix are the eigenvectors 
of the generalized Hartree-Fock-Bogoliubov (HFB) equa- 
tions (see ref. [f4|-^6[ for details). As such are, of course, 
determined only up to an overall phase (for each one). 

If the p-n pairing interaction is switched on, the angu- 
lar momentum-coupled phonon operator takes the form 
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Q?u = E™«< 6 ' J ) B t2(a, b, J, M) 

a, b 

+ Y$(a,b, J)B 12 (a,b,J,M)} 

+ j2i x Z( a > b > J ) B U a > b > J > M ) 



a<b 

fM—1,2 



+ Y™(a,b,J)B^(a,b,J,M)}, (57) 

The amplitudes X m , Y m are obtained by solving the 
QRPA eigenproblem 



A B\ X 



B A \ Y 



m I — J "m 



1 o \ f x r 

-1 I [ Y n 



with 





/A 11 ' 11 


411,22 


A ll,12 


A=\ 


A 22 ' 11 


,122,22 


^22,12 




\A 12 ' n 


^12,22 


^12,12 






£11,22 


£11,12 




£22,11 


£22,22 


£22,12 




KB 12 ' 11 


^12,22 


£12,12 



and 



X m = 




Y m = Kg 

Vk 



(58) 



(59) 



(60) 



(61) 



12 



The QRPA equation in Eq. ( pq ) represents a general 
equation for all excited states of a given even-even nu- 
cleus. The explicit form of the A and B submatrices 
has been reviewed in our previous publications |^4| ff"o| |. 
We note that if proton-neutron pairing interaction is ne- 
glected, the QRPA submatrices A 12 ' 11 , A 11 ' 12 , A 12 ' 22 
A 22 ' 12 and B 12 ' 11 , B 11 ' 12 , B 12 ' 22 , B 22A2 are equal to 
zero and the QRPA equation posesses two types of eigen- 
states: The eigenstates of the first type ( I ) with 
wave functions X m = (XJ},Xg,0), Y m = (Yff,Kj£,0) 
having the origin in proton-proton and neutron-neutron 
quasiparticle excitations. The eigenstate of the second 
type ( II ) with the wave functions X m = (0,0, Xg), 
Y m = (0, 0, YI2) generated by the phonon operator in 

Eqs. (|53|), and (|54]) are generalized as follows 



< mJ || u k (a+ a nb ) || 0+ >= <W2J+1 
x ^ rn{p,a,vb)[u lip {a)v vn {b)X^ l ,{a,b,J) 



/i,i/=l,2 



+ v w (a)u vn {b)Y™(a, b, J)], 



(62) 



< 0+ || u k (a+ a nd ) || mJ >= 4jV2J+T 
x ^ m {V a > u b)[v w {c)u vn {d)X™ v (c,d,J) 

(U,i/=1,2 

+ « w (cK n (<i)y;(c I d 1 j)], (63) 



with 

m(/za, = (1 + (-l) J <5 Al ^ a6 )/(l + S^S ab )^ 2 . (64) 

We note that the X^(a, b, J) and Y™(a, b, J) amplitudes 
are calculated by the QRPA equation in Eq. (|5|) only for 
the configurations fia < vb ( i.e., fi = v and the orbitals 
are ordered a < b and /Lt = 1, v = 2 and the orbitals 
are not ordered). For different configurations X™ v (a, b, J) 
and Y™(a, b, J) in Eqs. (|62|), and ( |63] ) are given as follows: 

X- (a, 6, J) = -(-1)^- J X- (6, a, J), (65) 

F-(a, 6, J) = -(-l)j*+Wy-(6, a, J). (66) 

Note that in the limit in which there is no proton-neutron 
pairing, i.e., u 2p = i'2v = u in = v ln = 0, Eqs. (|||), and 
(p3|) reduce to Eqs. (m), and (p3) respectively by setting 



Ulp 



u p ,vi p 



v p ,u 2n = u n , and v 2n 



Clearly, in 



the case without proton-neutron pairing the eigenstates 
of type I do not contribute to the beta decay transition 
matrix elements in Eqs. (|62[), and (|63|). 
The overlap integral becomes 

<mJ\mJ>= ]T [X™(a,b,J)Xf v (a,b,J) 

-Y£(a,b,J)Y%(a,b,J)]. (67) 



Y" 



the above 



By setting X™ = Yfi = X™ 
overlap is reduced to that of Eq. ( |55| ) . 

To complete the discussion we mention that the single 
particle wave functions and energies were obtained by 
using a Coulomb corrected Woods Saxon potential. The 
interaction employed was the Brueckner G-matrix which 
is a solution of the Bethe-Goldstone equation employed 
using the Bonn OBEP . Proton and neutron number 
conservation in the initial and final state was respected 
on the average with 



(iV n p- < A np >)/A np < 10" 



(68) 



The BCS pp and nn parameters d pp and d nn were ob- 
tained by fits to the experimental proton and neutron 
gaps as in ref. ]4^-|46]|. The np strength parameter d np 
is fixed by a renormalization of the T = 1 J = pairing 
force as in ref. Jl4]-[l6| . In the QRPA calculations it is nec- 
essary to introduce renormalization parameters g pp and 
g p h for the particle-particle and particle-hole interactions, 
which in principle should be close to unity. Our adopted 
values were <7 PP = 1.0 and g p h = 0.8. For higher value of 
g p h the particle-hole interaction for some multipolarities 
is too strong. The lowest eigenvalue becomes imaginary 
and leads to a collapse of the correlated ground state. By 
the method outlined above we obtained the matrix ele- 
ments M G U T , xf, xh, Xf, Xgt, Xt> Xf^> Xgtu, Xp and 
X R for the nuclei 48 Ca, 76 Ge, 82 Se, 96 Zr, 100 Mo, 116 Cd, 
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128 Te, 130 Te and 136 Xe for the orbitals shown in Table 
H These matrix elements are given in Table 0. For com- 
parison we also present in the same Table Othe values 
obtained without p-n pairing. These last results differ 
slightly from those of our earlier calculations J26|-^8| due 
to the different model space employed. 

By glancing at the Table || we see that the effect of 
the inclusion of the p-n pairing is significant. Perhaps 
the most important matrix element is |Mq^(1 — xf)|, 
which connect us directly with the effective neutrino mass 
| < m„ > |. We see that the inclusion of the p-n pair- 
ing reduces the value of this matrix element. The largest 
reductions of \M& T (1 - X f)\ by factors 30.2, 6.7, 3.4, 
2.8 and 2.3 are associated with A = 100, 96, 48, 128 
and 76 systems. By these factors also the lower limits 
on | < m„ > | are enhanced in respect to the calcula- 
tions without p-n pairing. To study the influence of p-n 
pairing on the evaluation of the limits on lepton num- 
ber non-conserving parameters of right-handed currents 
it is necessary to calculate the integrated kinematical fac- 
tors G ok m& The y are listed in Table H3- A sma11 

differences with the values of ref. jl9| has origin in dif- 
ferent adopted value of nuclear radius R [see Eq. (E8|)]. 
The problem of the extraction of the lepton number non- 
conserving parameters we shall study in Sec. IV. 

In the nuclear systems A — 96, 100 and 116 we have no- 
ticed a sensitivity of this matrix element with respect to 
the number of orbitals employed. Since including in our 
present calculation with p-n pairing all the 15 orbitals 
employed in the earlier calculations was prohibitive in 
term of computer time, we decided to employ 12 orbitals. 
This is admissible since we are interested in compair- 
ing the results with and without p-n pairing in the same 
model space. Furthermore our present results without p- 
n pairing for 100 Mo agree with those of ref. || which used 
the same model space. The results so obtained for this 
nucleus are comparable with those of the other nuclei, 
which makes our choice reasonable. The other matrix 
elements depend a bit more strongly on the p-n pairing 
correlations. The effect is even stronger for some indi- 
vidual multipoles, especially if the corresponding matrix 
element is suppressed. 



IV. DECAY RATES 

The Of/3/3-decay can be expressed in terms of the lep- 
ton violating parameter < m v > /m e etc. defined in Sec. 
II as follows |27), namely 

P^r 1 = Om^I 2 {|V L | 2 + |V R | 2 - C[X L X R 
+C* 2 |A|V L co S V;i + C 3 \ri\X h costl>2 + G 4 |A| 2 + C 5 \r,\ 2 
+C 6 \\\\r ] \cos(il> 1 - i> 2 ) + Re{C 2 XX R + C 3 r)X R )} , (69) 

where X\, and X R are defined in Eqs. Q), and (|35|). Vi 
and ip2 are the relative phases between Xl and A and 



Vl and 77 respectively. The ellipses {...} indicate con- 
tributions arising from other particles, e.g., intermediate 
SUSY particles |48| or unusual particles which are pre- 
dicted by superstring models jl9| or exotic Higgs scalars 
|0| etc. 

The quantities Gqi are calculated using the prescrip- 
tion of ref. Q. The coefficients C[,Ci,i = 2 — 6 are 
combinations of kinematical functions and the nuclear 
matrix elements discussed in the previous section. They 
arc defined as follows. 



C 2 

g 3 
G 5 = 

C 6 = 



-(1 - Xf)(x 2 - G 03 ~ Xi+G 04 ), 

— (1 - Xf)(X2+ < -^03 — X1-G04 — x'pGq5 + XrGog), 

1-2 
X2-G02 + gXi+G 4 - gXi+X2-G 03 , 



X2+C02 

-2[X2-X2+ G 02 



1 2 ~ 2 

gXl-G04 - gXl-X2+G03 



(Xp) 2 G'o8 



XpXrGW + XrGo9), 



g(Xi+X 2 + + X2-Xi-)G 3 

+ 9X1+X1-G04]. (70) 



Here C[ S 10(e§ + 6e + 6)/(e 4 + lOe 3 , + 10e§ + 60e + 30), 
eo is the available energy in electron mass units. C[ is 
less than 10% and it can be safely neglected. 
The quantities Goi are defined as follows: 



G*oi — Goi/Goi (i — 2, 3, 4), 
G05 = 2G05/G01, 

G06 = jWe-RoGoe/Goi, 

G07 = 2(-rn i?o)Go7/Goi, 

G08 = 4Go8/Goi, 

G09 = (— w c -Ro) 2 Go9/Goi. 



(71) 



The values of the parameters Go;,i = 2,..., 6 are pre- 
sented in Table IV. The coefficients Gi,i = 2, ...,6 with 
and without p-n pairing are shown in Table [v|. 
The most stringent experimental limits are: 



A = 


48 : 


rpOv 
J l/2 


> 9.5 x 10 21 y (76%GL) [5], 


A = 


76 : 


1 l/2 


> 5.6 x 10 24 y (90%GL) [2], 


A = 


82 : 


J l/2 


> 2.7 x 10 22 y (68%GL) [3], 


A = 


96 : 


rp&\\ 

1 l/2 


> 3.9 x 10 19 y (geochem.) [11], 


A = 


100 


1 l/2 


> 4.4 x 10 22 y (68%GL) [4], 


A = 


116 


1 l/2 


> 2.9 x 10 22 y (90%GL) [6], 


A = 


128 


J l/2 


> 7.7 x 10 24 y {geochem.) [7] 


A = 


130 


1 l/2 


> 2.3 x 10 22 y (68%GL) [8], 


A = 


136 


1 l/2 


> 3.4 x 10 23 y (68%GL) [9]. 
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Using Eq. fl69|), the functions Goi of Table [V and the 
matrix elements Mq^ of Table || with p-n pairing we get 
the constrains: 



A — 




/ 

V 


1 < 9K v 10" 

.j _^ O.U A _LU 


-9 

3 


4 — 


I . 


r 

V 


1 <r r v 1 n~ 

.J- _\ u.D X 1U 


12 

3 


4 — 


oz . 


r 

1- 


1 / 7Q Y 10" 


10 

3 


A = 


96 : 


{• 


.} < 4.4 x 10" 


-6 

3 


A = 


100 


{■ 


..} < 1.2 x 10" 


-9 

3 


A = 


116 


{■ 


..} < 3.9 x 10" 


-8 

3 


A = 


128 


{■ 


..} < 3.8 x 10" 


-11 


A = 


130 


{• 


..} < 2.3 x 10" 


-10 


A = 


136 


{• 


..} < 1.4 x 10" 


-11 



Here {...} indicates the quantity which is enclosed in the 
curly braket of Eq. j69|) . Clearly the nucleus with the 
smallest value of {...} is going to provide the most strin- 
gent limit on the lepton violating parameters. The life 
time itself is not a clear indicator since the function Goi 
varies from nucleus to nucleus. Large Goi, i.e., large 
phase space, leads to short life times for a given lepton 
violation parameter. Thus the most stringent limits are 
expected from the A=76 system. 

To impose limits on Xl, -Xr, A, r\ one must make four - 
dimensional plots making some assumptions about ipi , ip2 
and the relative signs of A and rj with \R- Then for a 
given value of Xl one can extract limits on < m„ > 
+ m o < Sn >l X H /(X F - 1 )- To extract limits on 
< m v > and < -p— > one must make further plots 
(knowing Xh> Xf horn the calculations). This is re- 
ally a complicated procedure to be worth doing only 
if and when Oi//3/3-decay is definitely seen. At present 
we will constrain the above parameters by assuming 
that one mechanism at a time dominates. The limits 
thus obtained appear in Table [V]|. We must mention 
that in the case of heavy neutrino only the parameters 
« jfc > R (e 2 + « 2 ))- 1 =< 

this way. The parameter < 



l 

M K 



1 

/2 



> 



> L can be extracted 
-l 

R 

2 _ 



shown in Table VI 

was obtained by taking |L8|] e 2 + n 2 = 10~ 2 . In line with 
what we mentioned above the extraction of the parame- 
ter r) depends on G5 alone. In all cases G5 is dominated 
by Xp and /or Xr ■ So in Table VII we present two values 
of 77, one with nuclear recoil included and one without 
recoil. With the possible exception of the A = 100 and 
A = 128, 130 system, the recoil contribution is dominant. 
In the Te isotope x'p and Xr compete with each other. 

Another lepton violating process is the Qv/30- decay 
with Majoron emission. The corresponding expression is 



XF) 



(72) 



where 



r, _ 1 2 ffiM 

G °*°-(2.) 2£ %( eo r 



5o (eo) = 4 + W4 + 40e 2 + 60e + 30, 
5 i (eo) = 4 + U4 + 84e 2 + 210e + 210, 
1 



Err( 11 )rr( 11 ) 



V2 



e i(ai - Qi) gij. 



(73) 



e is the available energy in units of the electron mass. 
<7ij is the coupling of the Majoron to the neutrino mass 
eigenstates, i.e., 



£ = ^P iLl / jRX ° + h.c. 



It can also be written as 
9, 



r 9cc c 



h.c. 



The corresponding experimental limits are 

" > 7.2 x 10 20 (90%GL) [14], 

^Jf > 3.9 x 10 22 ( 90%GL) [2], 

. l > 1.6 x 10 21 (68%GL) [3], 



(74) 



(75) 



A = 


48 : 


ml 
1 1 


A = 


76 : 


1 1 


A = 


82 : 


rpl 
1 1 


A = 


96 : 


rpl 
1 1 


A = 


100 


T[ 


A = 


116 


T[ 


A = 


128 


T[ 


A = 


130 


T{ 


A = 


136 : 


t 



/2 



,19 



100 : T"J 2 X > 7.9 x 10 2U (68%GL) [15], 
°^ x ° > 1.8 x 10 19 (99%GL) [6], 

all ^ -7 7 w in24 /„„„>„„, \ M 



all 



0v,x u 



From the above experimental limits and the values of 
Goi, x o (Table and |Mg^.(l-xp)| (Tableg) we obtain 
the limits of |ry x o| listed below: 



.4 
A 
A 
A 



76 
82 
96 



A = 100 
A = 116 
A = 128 
A = 130 
A = 136 



|?7 x o| < 7.0 x 10~ 4 , 
| % o| < 1.4 x 10- 4 , 
|t7 x o| < 1.9 x 10~ 4 , 
|t7 x o| < 3.6 x 10~ 3 , 

< 9.9 x 10~ 3 , 

< 2.6 x 10~ 3 , 



\V X ° 
\V X ° 

K 



< 5.7 x 10"°, 

< 1.5 x 10" 4 , 

< 1.3 x 10" 4 . 



V. SUMMARY AND CONCLUSIONS 

In the present work we have evaluated the nuclear ma- 
trix elements entering the double beta decay of the ex- 
perimentally most interesting nuclear systems. We have 
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employed the Quasiparticle Random Phase Approxima- 
tion which seems to be the most practical method of nu- 
clear structure calculation of nuclear systems which are 
far away from closed shells. In these calculations we have 
included the proton-neutron pairing correlations which 
have been neglected in the previous calculations. Wc 
have found that such correlations have important effects 
on all the needed matrix elements and should not be ne- 
glected. The magnitude of the effect depends, of course, 
on the type of operator employed, i.e., on the mechanism 
for the 0i//3/?-decay. We will concentrate our discussion 
on those matrix elements which were not unusually sup- 
pressed. We will give their value both without p-n pair- 
ing correlations and after such correlations are turned on. 
Since we have assumed that one mechanism at a time is 
important we can summarize our results as follows. 

i) Light neutrino mass mechanism. The relevant nu- 
clear matrix element is \Mq^(1 — x F )| 2 . The five largest 
matrix elements are 9.2, 6.2, 5.8, 5.6 and 4.9 associated 
with A = 76, 128, 96, 130 and 82 respectively. Once p-n 
pairing correlations are turned on they become 1.8, 0.77, 
0.13, 2.2 and 2.7 respectively, i.e., we have a reduction 
factor ranging between 2 and 8. In the most interesting 
case of the A — 76 system we have a reduction of about 
5. The effect on the extraction on the neutrino mass is 
less pronounced. (The square root of the above factor). 

ii) Heavy intermediate neutrino mass mechanism. The 
relevant nuclear matrix element is |x h Mq^,| 2 . The five 
largest matrix elements are 4.0 x 10 4 , 1.1 x 10 4 , 1.0 x 10 4 , 
0.98 x 10 4 and 0.96 x 10 4 for A = 76,48, 128,96 and 82 
respectively. With p-n pairing correlations they become 
0.37 x 10 4 , 6.0, 0.19 x 10 4 , 1.4 x 10 2 and 0.43 x 10 4 . Notice 
the almost complete suppression for A = 48 and 96 and 
the large reduction factor for A — 76 (about 10). 

iii) The mass independent A- mechanism. The rele- 
vant matrix element is C 4 |Mq^| 2 . The largest matrix el- 
ements 14.5, 14.1, 13.6, 13.4 and 6.2 are associated with 
A = 96, 76, 82, 130 and 136 respectively. They become 
0.49, 3.07, 6.5, 8.1 and 3.8 respectively. Notice that quite 
unexpectedly the matrix element of the A = 128 system 
is much smaller than that for the A — 130 system. In the 
case of the A = 76 we have a reduction factor of about 
5. 

iv) The neutrino mass independent r\- mechanism. The 
relevant matrix element is C 5 | | 2 . The largest values 
are 4.1 x 10 5 , 3.7 x 10 5 , 2.9 x 10 5 , 1.4 x 10 5 and 9.1 x 10 4 
associated with A — 76, 128, 130, 136 and 96. They are 
reduced to 5.5 x 10 4 , 2.4 x 10 5 , 1.7 x 10 5 , 1.8 x 10 5 and 
1.7 x 10 4 . In the case of the A = 76 system the reduction 
factor is about 7. 

We do not fully understand why the effect of the p-n 
pairing correlations should be so large. In the case of 
A = 100 this effect is much more dramatic. On the con- 
trary for the A = 136 system the matrix elements are 
fairly large and not much effected by such correlations. 
Finally we do not understand why the effect of p-n cor- 
relations is so different on the two Tellurium isotopes. 



From Eq. (J69|) it is clear that the A = 76 system pro- 
vides the most stringent limits on the lepton violating 
parameters. This is partly due to the large matrix el- 
ements obtained for this system but mainly due to the 
fact that the experimental life-time limit is the best. Un- 
fortunately, the introduction of p-n pairing correlations 
makes the extracted limits of the lepton violating param- 
eters less stringent. In fact we find 



< m„ > | < 



eV 
eV 



< 



1 
1 



>T 



> 



1-1 



> 



< A > | < 

<n > I < 



0.8 
1.8 

3.8 x 10 7 GeV 
1.6 x 10 7 GeV 

3.8 x 10 5 GeV 

1.6 x 10 5 GeV 

1.3 x 10~ 6 

2.7 x 10~ 6 ' 

7.4 x 10" 9 
2.0 x 10~ 8 ' 

8.4 x 10" 5 

1.9 x 10~ 4 " 



In the above expressions the upper (lower) values corre- 
spond to the case without (with) p-n pairing correlations. 
It is interesting to note that the A — 128 system, in spite 
of the fact that its 0v(3f3 the decay width is kinemati- 
cally suppressed, provides quite stringent limits on the 
lepton violating parameters with the possible exception 
of A. This is quite surprising since the nuclear matrix el- 
ements involved are not favored compared to the A = 76 
system and in any case they should not be very differ- 
ent from those of the A — 130 system. Furthermore the 
extracted limits on the lepton violating parameters will 
become even more stringent if the 0v life time is used 
since Ty 2 > T^j l 2 . We are, therefore, inclined to suspect 

that the life time of T®" 2 > 5.6 x 10 24 years is quite a bit 
exaggerated. We will not, however, elaborate further on 
this controversial point concerning the life time of this 
long lived isotope. 

Finally we should mention that the above extracted 
limits still suffer from uncertainties of nuclear origin. We 
should not forget that the effect of the interference be- 
tween the various mechanisms has not been taken into 
account. 
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TABLE I. The single particle orbitals used in the present calculations. 



nuclear 48 Ca 76 Ge 82 Se 96 Zr 100 Mo llfi Cd 128 Te 130 Te ui& Xe 

orbitals 



1 


OSi/2 


lSl/2 


lSl/2 


0/7/2 


0/7/2 


0/5/2 


1P3/2 


1P3/2 


1P3/2 


2 


0^3/2 


0/7/2 


0/7/2 


0/5/2 


0/5/2 


1P3/2 


lPl/2 


lPl/2 


lPl/2 


3 


0pi/2 


0/5/2 


0/5/2 


1P3/2 


1P3/2 


lPl/2 


057/2 


0ff7/2 


0ff7/ 2 


4 


0d 5 / 2 


lf>3/2 


1P3/2 


lPl/2 


lPl/2 


059/2 


0ff7/2 


0ff7/2 


0ff7/ 2 


5 


0d 3 / 2 


lPl/2 


lPl/2 


0fl9/2 


0fl9/2 


0ff7/ 2 


14/2 


14/2 


14/2 


6 


lSl/2 


059/2 


0ff9/2 


0ff7/ 2 


0ff7/2 


14/2 


14/2 


W3/2 


14/2 


7 


0/7/2 


057/2 


057/2 


14/2 


14/2 


1^3/2 


2si/ 2 


2S1/2 


2S1/2 


8 


0/5/2 


14/2 


14/2 


14/2 


14/2 


2«l/2 


0/lll/2 


0^11/2 


Oft-11/2 


9 


1P3/2 


1^3/2 


14/2 


2S1/2 


2S1/2 


0^11/2 


04/2 


04/2 


O/lg/2 


10 


lPl/2 


2S1/2 


2Si/ 2 


0^11/2 


0^11/2 


0h 9 /2 


I/7/2 


1/7/2 


1/7/2 


11 








O/lg/2 


04/2 


1/7/2 








12 












1/5/2 
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TABLE II. The matrix elements of Oi//3/3-decay for 48 Ca, 76 Ge, 82 Se, 96 Zr, 100 Mo, 116 Cd, 128 Te, 130 Te and 136 Xe calculated 
in the framework of QRPA with and without p-n pairing. 



nucleus 


48 Ca 


76 Ge 


82 5e 


96 Zr 


100 Mo 


116 Cd 


128 Te 


130 Te 


136 Ae 








QRPA without p-n pairing 










Mgt 
Xf 

\Mgt(1-Xf)\ 
Xh 


-0.785 
-0.468 
1.152 
-134.9 


2.929 
-0.038 
3.040 
-68.37 


-2.212 
-0.008 
2.230 
-44.27 


2.097 
-0.149 
2.409 
-47.24 


0.615 
-0.766 
1.086 
-124.8 


0.449 
-1.103 
0.944 
-47.06 


2.437 
-0.0179 
2.480 
-41.54 


2.327 
-0.004 
2.335 
-39.82 


1.598 
0.028 
1.553 
-21.92 


Xf 

ACj 1 

Xt 

XFu 

Xgtw 
Xp 


-0.504 
0.975 
-0.212 
-0.437 
1.057 
0.168 


-0.035 
1.077 
0.244 
-0.038 
0.916 
-1.147 


-0.004 
1.050 
0.079 
-0.013 
0.960 
-0.049 


-0.168 
1.143 
0.121 
-0.130 
0.845 
-0.836 


-0.817 
1.174 
-0.477 
-0.709 
0.683 
-3.843 


-1.173 
1.074 
-0.812 
-1.032 
0.859 
-3.891 


-0.022 
1.097 
0.307 
-0.012 
0.894 
-1.400 


-0.007 
1.097 
0.282 
0.001 
0.895 
-1.451 


0.022 
1.123 
0.349 
0.036 
0.875 
-1.627 


Xr 


172.1 


193.0 


124.2 


113.8 


105.1 


-151.5 


157.1 


149.0 


124.8 








QRPA with p-n 


pairing 










XF 

\Mqt(1-Xf)\ 

XH 


-0.405 
0.158 
0.341 
6.075 


1.846 
0.274 
1.340 
-32.75 


-1.153 
-0.416 
1.633 
-57.20 


0.280 
2.282 
0.358 
-41.64 


-0.584 
0.939 
0.036 
-14.22 


0.119 
-6.784 
0.926 
-453.8 


1.270 
0.308 
0.879 
-34.02 


1.833 
0.184 
1.495 
-55.72 


1.346 
0.066 
1.257 
-35.37 


XF 
XGT 

XT 
XF^ 

Xgtlj 
Xp 


0.184 

1 99fi 

0.130 
0.131 
0.775 
-0.009 


0.322 
1 1 94 

0.214 
0.235 
0.876 
-0.479 


-0.467 
1.082 
0.179 
-0.379 
0.927 
-1.621 


2.601 
1.587 
0.209 
2.069 
0.335 
-4.802 


1.067 
0.934 
0.853 
0.812 
1.142 
2.519 


-7.400 

n Q97 

-3.991 
-6.170 
0.938 
-7.592 


0.370 

1 1 t\Q 

0.343 
0.260 
0.831 
-2.907 


0.218 
i lie; 

0.411 
0.159 
0.879 
-0.993 


0.082 

1 1 R7 

0.332 
0.052 
0.832 
-2.441 


XR 


57.32 


129.3 


131.1 


157.3 


162.2 


-333.6 


158.6 


192.6 


138.4 
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TABLE III. The integrated kinematical factors Gok for + — > + transition of {(3f3)o v — decay. The definition of Got is given 
in ref. [19] in Eqs. (3.5.17-21). 



(/3/3)ov — decay : + — > + transition 







48 Ca 


7e Ge 


82 Se 


96 Zr 


100 Mo 


116 Gd 


128 Te 




130 Te 


136 Xe 


(Ei- 


- Ef) [MeV] 


5.294 


3.067 


4.027 


4.372 


4.055 


3.830 


1.891 




3.555 


3.503 


Goi 


[io^V 1 ] 


8.031 


0.7928 


3.524 


7.362 


5.731 


6.233 


2.207x10" 


-l 


5.543 


5.914 


G02 


[io^V 1 ] 


5.235 


0.1296 


1.221 


3.173 


2.056 


1.957 


6.309x10" 


■3 


1.441 


1.483 


G(i3 


[W-^y- 1 ] 


6.037 


0.4376 


2.413 


5.380 


4.036 


4.305 


6.177x10" 


-2 


3.669 


3.890 


Go4 


[io^V 1 ] 


1.705 


0.1538 


0.724 


1.530 


1.178 


1.269 


3.368x10" 


-2 


1.113 


1.183 


Gor, 


[io^V 1 ] 


1.265 


0.253 


0.931 


2.009 


1.718 


2.118 


1.390x10" 


-1 


2.083 


2.298 


Goe 


[io^V 1 ] 


1.398 


0.196 


0.665 


1.226 


1.009 


1.103 


6.969x10" 


-2 


1.011 


1.077 


Go7 


[lO^V 1 ] 


11.46 


1.248 


5.523 


12.07 


9.563 


10.69 


4.363x10" 


-1 


9.544 


10.25 


Gos 


[io^V 1 ] 


5.247 


0.793 


3.852 


9.886 


8.109 


10.20 


4.227x10" 


-1 


9.749 


10.84 


G09 


[W-Vy- 1 ] 


6.262 


0.491 


1.980 


3.686 


2.819 


2.800 


1.125x10" 


-1 


2.335 


2.424 



Goi,x° [lO" 14 */- 1 ] 


2.425 


0.0763 


0.6202 1.5315 


1.0230 


0.9879 


5.206xl0~ 3 


0.7487 


0.7734 


TABLE IV. The kinematical functions Goi, i = 2 - 9. 


They are § 


,dven in the notation of Pantis et al. [27] . 




nuclear 


G02 


Go3 


Go4 


Gl)5 


Goe 


Go 7 


Gl)8 


G09 


transition 


















i* C a - i8 Ti 


6.518 


0.752 


0.212 


31.50 


0.450 


73.87 


2613. 


0.522 




1.635 


0.552 


0.194 


63.93 


0.745 


95.01 


4001. 


0.563 


82 Se -» 82 Kr 


3.465 


0.685 


0.205 


52.9 


0.584 


96.99 


4372. 


0.538 


96 Zr 96 Mo 


4.310 


0.731 


0.208 


54.58 


0.543 


107.0 


5371. 


0.5324 


ioo Mo _ ioo Ru 


3.588 


0.704 


0.206 


59.97 


0.582 


110.3 


5660. 


0.5375 




3.140 


0.691 


0.204 


67.94 


0.614 


119.1 


6547. 


0.5419 


128 Te ^ 128 Xe 


0.286 


0.280 


0.153 


126.0 


1.133 


141.9 


7662. 


0.6565 


130 Tg ^ 130 Xe 


2.599 


0.662 


0.201 


75.15 


0.658 


124.2 


7035. 


0.5483 


136 Xe _ 136 Ba 


2.507 


0.658 


0.200 


77.72 


0.667 


127.0 


7331. 


0.5497 
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TABLE V. The coefficients C\, i = 1, 2, 3, 4, 5, 6 which are combinations of nuclear matrix elements and phase space factors 



[see Eq. (55) of the text] needed for the extraction of | < m v > |, A, r\ etc. from the data. 

















nuclear 


Gi 


g 2 


c 3 




c 5 


g 6 


transition 


















QRPA without p-n pairing 








^Ca -> 48 Ti 


2.154 


-0.96 


-1.05 x 10 2 


7.43 


1.34 x 10 4 


-6.92 




1.078 


-0.66 


-2.26 x 10 2 


1.65 


4.73 x 10 4 


-3.10 


82 Se -> 82 Kr 


1.015 


-0.51 


-7.61 x 10 


2.78 


8 90 x 1 3 


-5.42 


96 Zr 96 Mo 


1.321 


-0.75 


-1.23 x 10 2 


3.29 


2.08 x 10 4 


-5.32 


ioo Mo ^ 100^ 


3.118 


-0.26 


-5.12 x 10 2 


1.51 


1.34 x 10 5 


1.13 


™Cd -> 116 


4.421 


-0.18 


-3.55 x 10 2 


2.03 


4.13 x 10 4 


2.35 


128 Te ^ 128 Xe 


1.036 


-0.40 


-3.61 x 10 2 


0.33 


6.24 x 10 4 


-0.65 


130 Te ^ 130 Xe 


1.007 


-0.76 


-2.08 x 10 2 


2.48 


5.38 x 10 4 


-4.98 


136 Xf , ^ 136 Ba 


0.944 


-0.78 


-2.04 x 10 2 


2.43 


5.37 x 10 4 


-5.16 








QRPA with p-n pairing 










0.708 


-0.24 


-2.25 x 10 


2.72 


1.76 x 10 3 


-6.67 


76 Ge -» 76 Se 


0.526 


-0.19 


-9.27 x 10 


0.90 


1.62 x 10 4 


-2.53 


82 Se _ 82^ 


2.005 


-1.52 


-2.30 x 10 2 


4.86 


4.13 x 10 4 


-5.41 


96 Zr -> 96 Mo 


1.643 


-3.01 


4.48 x 10 2 


6.24 


2.18 x 10 5 


17.6 


ioo Mo _ ioo fiw 


0.0037 


-0.062 


3.27 


4.95 


4.99 x 10 3 


-19.0 


116 Cd 116 Sra 


60.59 


-5.86 


-2.35 x 10 3 


10.4 


1.36 x 10 5 


70.8 


128 Te > 128 Xe 


0.479 


-0.13 


-3.78 x 10 2 


0.17 


1.47 x 10 5 


-0.56 


130 Te ^ 130 Xe 


0.665 


-0.62 


-1.65 x 10 2 


2.41 


5.11 x 10 4 


-5.91 


136 Xe ^ 136 Ba 


0.872 


-0.66 


-2.64 x 10 2 


2.11 


9.71 x 10 4 


-4.53 



TABLE VI. The limits on lepton number non-conserving parameters | < ra„ > j, | < M 1 > | 1 , and | < \ 1 deduced 

from the experimental limits of Quj3f3 decay lifetimes for the nuclei studied in this work. For the extraction of the parameter 
| < -p— > H [ the value of e 2 + k 2 = 10 -2 was assumed. Only one term was assumed dominant. T®J^~ cV is the calculated 
Qvj3j3 half-life times assuming | < m v > | = leV. 

nucleus \<m v >\ \ < -fa > L p 1 | < > R p 1 T%?*» 

[eV] [GeV] [GeV] [years] [years ] ref. 



QRPA without p-n pairing 



AS Ca 


< 16 


> 1.9 x 10 6 


> 1.9 x 10 4 


2.5 x 10 24 


76 Ge 


< 0.8 


> 3.8 x 10 7 


> 3.8 x 10 5 


3.6 x 10 24 


S2 Se 


< 7.4 


> 2.8 x 10 6 


> 2.8 x 10 4 


1.5 x 10 24 


96 Zr 


< 125 


> 1.4 x 10 5 


> 1.4 x 10 3 


6.1 x 10 23 


100 Mo 


< 9.4 


> 2.0 x 10 6 


> 2.0 x 10 4 


3.9 x 10 24 


116 Cd 


< 13 


> 4.0 x 10 5 


> 4.0 x 10 3 


4.7 x 10 24 


128 Te 


< 1.6 


> 1.2 x 10 7 


> 1.2 x 10 5 


1.9 x 10 25 


130 Te 


< 6.1 


> 3.1 x 10 6 


> 3.1 x 10 4 


8.6 x 10 23 


136 Xe 


< 1-7 


> 6.6 x 10 6 


> 6.6 x 10 4 


1.8 x 10 24 



> 9.5 x 10 21 

> 5.6 x 10 24 

> 2.7 x 10 22 

> 3.9 x 10 19 j 

> 4.4 x 10 22 

> 2.9 x 10 22 

> 7.3 x 10 24 

> 2.3 x 10 22 j 

> 6.4 x 10 23 



QRPA with p-n pairing 



48 Ca 


< 54 


> 7.6 x 10 4 


> 7.6 x 10 2 


2.8 x 10 25 


76 Ge 


< 1.8 


> 1.6 x 10 7 


> 1.6 x 10 5 


1.8 x 10 25 


S2 Se 


< 10 


> 1.3 x 10 6 


> 1.3 x 10 4 


2.8 x 10 24 


96 Zr 


< 841 


> 1.4 x 10 4 


> 1.4 x 10 2 


2.7 x 10 25 


100 Mo 


< 285 


> 6.4 x 10 6 


> 6.4 x 10 4 


3.6 x 10 27 


116 Cd 


< 13 


> 2.8 x 10 5 


> 2.8 x 10 3 


4.9 x 10 24 


128 Te 


< 4.6 


> 7.4 x 10 6 


> 7.4 x 10 4 


1.5 x 10 26 


130 Te 


< 9.6 


> 4.2 x 10 6 


> 4.2 x 10 4 


2.1 x 10 24 


136 Xe 


< 2.1 


> 9.3 x 10 6 


> 9.3 x 10 4 


2.8 x 10 24 



10 



> 9.5 x 10 2 

> 5.6 x 10 2 

> 2.7 x 10 2 

> 3.9 x 10 

> 4.4 x 10 2 

> 2.9 x 10 2 

> 7.3 x 10 2 

> 2.3 x 10 25 

> 6.4 x 10 2 



s 



15 



TABLE VII. The limits on lepton number non-conserving parameters | < A > and | < f] > | deduced from the experimental 
limits of Ov/3/3 decay lifetimes for the nuclei studied in this work. For < rj > our results with and without inclusion of the recoil 
nuclear matrix elements are presented. Only one term was assumed dominant. 



nucleus | < A > | | < rj > 

without with 
recoil recoil 



T°J 2 e * p (Qu) 
[years] ref. 



QRPA without p-n pairing 



48 Ca 


< 1.7 x 1(T 5 


< 5.3 x 1(T 6 


< 4.0 x 10" 7 


76 Ge 


< 1.3 x 1(T 6 


< 2.2 x 1(T 8 


< 7.4 x 1(T 9 


82 Se 


< 8.8 x 1CT 6 


< 4.1 x 1(T 6 


< 1.5 x 1(T 7 


96 Zr 


< 1.6 x 1(T 4 


< 4.6 x 10~ 6 


< 1.9 x 1(T 6 


100 Mo 


< 2.6 x 1(T 5 


< 1.1 x 10~ 7 


< 8.8 x 10 -8 


116 Gd 


< 3.7 x 1(T 5 


< 1.7 x 10~ 7 


< 2.6 x 10~ 7 


128 Te 


< 5.6 x 1(T 6 


< 2.6 x 1(T 8 


< 1.3 x 1(T 8 


130 Te 


< 7.6 x 10~ 6 


< 9.9 x 1(T 8 


< 5.2 x 1(T 8 


136 Xe 


< 2.1 x 10~ 6 


< 2.3 x 10 -8 


< 1.4 x 1(T 8 



> 9.5 x 10 21 

> 5.6 x 10 24 

> 2.7 x 10 22 

> 3.9 x 10 19 

> 4.4 x 10 22 

> 2.9 x 10 22 

> 7.3 x 10 24 

> 2.3 x 10 22 

> 6.4 x 10 23 



QRPA with p-n pairing 



AS Ca 


< 5.4 x 10" 6 


< 4.3 x 1(T 5 


< 2.1 x 10~ 6 


76 Ge 


< 2.7 x 1(T 6 


< 8.5 x 1(T 8 


< 2.0 x 10~ 8 


82 5e 


< 1.3 x 1(T 5 


< 2.6 x 1(T 7 


< 1.4 x 10 -7 


96 Zr 


< 8.4 x 1(T 4 


< 6.0 x 10 -6 


< 4.5 x 10~ 6 


100 Mo 


< 1.5 x 10 -5 


< 1.8 x 10 -7 


< 4.8 x 10~ 7 


116 Gd 


< 6.1 x 1(T 5 


< 3.2 x 10~ 7 


< 5.4 x 10 -7 


128 Te 


< 1.5 x 1(T 5 


< 2.4 x 10~ 8 


< 1.6 x 10~ 8 


130 Te 


< 9.8 x 10~ 6 


< 1.8 x 10~ 7 


< 6.8 x 10~ 8 


136 Xe 


< 2.6 x 10 -6 


< 1.8 x 10 -8 


< 1.2 x 10~ 8 



111 



> 9.5 x 10 2 

> 5.6 x 10 2 

> 2.7 x 10 2 

> 3.9 x 10 

> 4.4 x 10 2 

> 2.9 x 10 2 

> 7.3 x 10 2 

> 2.3 x 10 25 

> 6.4 x 10 2 



1G 



